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INTRODUCTION 


As the era of the Space Transportation System progresses, and the possibility of 
a permanent manned presence in space unfolds , the necessity for improved radiation 
protection methods , for career astronauts who will spend appreciable periods of time 
in the space environment, increases in significance (ref, 1). The combination of 
long-duration missions and career astronauts will cause the high-energy heavy-ion 
(HZE) component of galactic cosmic rays to assume major radiobiological significance, 
especially for cumulative damage to nonregenerative tissues (ref. 2) and latent 
carcinogenesis (ref. 3). The greater ranges of these particles, when compared with 
typical thicknesses of spacecraft walls, suggest that nuclear attenuation and 
fragmentation may provide the only viable means of protection. Unfortunately, the 
secondary radiations produced by these galactic cosmic rays as they fragment also 
include HZE particles, so that appreciable thicknesses of bulk material may not sub- 
stantially reduce the transmitted doses (ref. 3). As a result, HZE particle shield- 
ing may prove to be impractical for space applications. The need remains, however, 
to develop an accurate fragmentation theory in order to adequately assess the effects 
of fragmentation by existing spacecraft structure and human body mass on the doses 
received by critical organs of interest. Alternatives to passive shielding methods, 
such as using electromagnetic fields as HZE particle shields, do not appear to be 
technologically (ref. 4) or economically feasible (ref. 5). However, no detailed 
study of their possible application to the entire spectrum of space radiations has 
yet been performed. 

In this work, a simple theory of HZE fragmentation based upon a two-step 
abrasion-ablation model is presented. For the abrasion step, a quantum-mechanical 
formalism based upon an optical model potential approximation to the exact nucleus- 
nucleus multiple-scattering series is used (refs. 6 to 9 ) . The charge dispersions 
of the excited projectile prefragments (the pieces of projectile nuclei remaining 
immediately after the abrasions) are calculated using both a hypergeometric distribu- 
tion and a method based upon the zero-point oscillations of the giant dipole reso- 
nance (ref. 10). The prefragment excitation energies are estimated from the geomet- 
ric "clean-cut" abrasion-ablation model of Bowman, Swiatecki, and Tsang (ref. 11). 

The ablation part of the fragmentation process is analyzed by calculating compound 
nucleus decay probabilities for the various particle emission channels using the 
Monte Carlo program EVAP-4, developed at Oak Ridge National Laboratory (ref. 12). 
Finally, to illustrate the theory, theoretical predictions of elemental production 
cross sections for 1 .88-GeV/nucleon iron projectiles colliding with stationary tar- 
get nuclei of carbon, silver, and lead are made and compared with experimental data 
(ref. 13) and with the predictions from the semiempirical relations of Silberberg 
et al. (ref. 14) . 


FRAGMENTATION THEORY 

In the abrasion-ablation fragmentation model, the projectile nuclei, moving at 
relativistic speeds, collide with stationary target nuclei. In the abrasion step, 
those portions of the nuclear volumes which overlap are sheared away by the colli- 
sion. The remaining projectile piece, called a prefragment or primary residue, con- 
tinues its trajectory with essentially its precollision velocity. As a result of the 
dynamics of the abrasion process, the prefragment is highly excited and subsequently 


decays by the emission of gamma radiation and/or nuclear particles. This step is the 
ablation stage. The resultant isotope, sometimes referred to as a secondary product, 
is the nuclear fragment whose cross section is measured. The abrasion process can be 
analyzed using classical geometric arguments (refs. 11 and 15) or methods obtained 
from formal quantum scattering theory (refs. 6 to 9 and 16). The ablation stage can 
be analyzed from geometric arguments (ref. 11) or more sophisticated methods based 
upon Monte Carlo or intranuclear cascade techniques (refs. 10, 12, 15, and 16). 
Predictions of fragmentation cross sections can also be made using the approximate, 
semiempirical, parameterization formulas of Silberberg et al. (ref. 14). 


Abrasion Cross Sections 

From reference 9, the cross section for abrading m projectile nucleons is 

a = ( A p) / d 2 S[l - P(b)f P(b) Ap (1) 

m ^m J J L 

where P(&) is the probability, as a function of impact parameter, for not removing 
a nucleon in the collision and Ap, the residual mass (prefragment) number, is 

A = A - m (2) 

F P 


Within the context of eikonal scattering theory, an optical model potential approxi- 
mation to the exact nucleus -nucleus multiple-scattering series, which includes Pauli 
correlation effects, yields (ref. 9) 


P(&) = exp[-A T a(e)l(S) ] 


(3) 


with 


1(b) = [ 2nB(e ) ] -3/2 / dz o / d 3 ! T p T (! T ) / d 3 y P p <S+* 0 +y+£ T > 
x exp[-y 2 /2B(e) ][l - C(y)] 


(4) 


(A list of symbols and abbreviations used in this paper appears after the refer- 
ences.) Methods for determining the appropriate nuclear distributions (i - P,T) 

and constituent-averaged nucleon-nucleon cross sections a(e) are described in 
reference 17. Values for the nucleon-nucleon scattering slope parameter are obtained 
from the parameterization, appropriate for diffractive scattering, given in refer- 
ence 18. From reference 9, the Pauli correlation is 


C(y) = 0.25 exp(-k 2 y 2 /10) 
K 


(5) 
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where 


k F = 1 ,36 fm -1 


Prefragment Charge Distributions 

Since the abraded nucleons consist of protons and neutrons, which are not iden- 
tical, a prescription for calculating the charge dispersions of the prefragments is 
needed in order to calculate final, isotope, and/or elemental production cross sec- 
tions caused by the fragmentation process. Two such methods are used in the fragmen- 
tation theory described in this work. The first method (ref. 19) treats the neutron 
and proton distributions as completely uncorrelated. The cross section for forming a 
particular prefragment of mass Aj and charge Zj is then given in terms of the 
cross section for abrading m nucleons a m as 


•wvv 



( 6 ) 


where z out of the original Z projectile nucleus protons are abraded along with 
n out of the original N projectile neutrons. Note that 


A = N + Z 
P 


(7) 


and 


m - n + z 


( 8 ) 


with 


Z. = Z - z 
J 


( 9 ) 


and 


A = Ap — m (10) 

The hypergeometric distribution is based on the assumption that there is no correla- 
tion at all between neutron and proton distributions. Therefore, unphysical results 
such as abrading all neutrons or protons from a nucleus while leaving the remaining 
fragment intact could occur. 
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As an alternative to the hyper geometric distribution , Morrissey et al. (ref. 10) 
proposed a charge dispersion model based upon the zero-point vibrations of the giant 
dipole resonance of the projectile nucleus. In this model, equation (6) becomes 


a , (Z . ,A. ) 

abr j j 



where the variance (dispersion) is 


a z = 2.61 9 (u/Ap ) 1 ( Z/Ap ) (dm/db) ( 1 + u) 3/ ^ 4 


(11 ) 


( 12 ) 


with 


u 


3J 


Q<V 


1/3 


(13) 


In the droplet model of the nucleus, the coefficients J and Q have the nominal 
values of 25.76 MeV and 11.9 MeV, respectively (ref. 10). The rate of change of the 
number of nucleons removed as a function of impact parameter (dm/db) is calculated 
numerically using the geometric abrasion model of reference 10. The normalization 
factor Nj insures that for a given value of Aj / the discrete sum over all allowed 
values of Zj yields unity for the dispersion probabilities. This overall normal- 
ization is a new feature of this work and is not included in the original model of 
reference 10. 


Prefragment Excitation Energies 

The excitation energy of the projectile prefragment following abrasion of m 
nucleons is calculated from the clean-cut abrasion formalism of references 11 and 15. 
For this model, the colliding nuclei are assumed to be uniform spheres of radii 

(i = P,T) • In the collision, the overlapping volumes shear off so that the resul- 
tant projectile prefragment is a sphere with a cylindrical hole gouged out of it. 

The excitation energy is then determined by calculating the difference in surface 
area between the misshapen sphere and a perfect sphere of equal volume. This excess 
surface area A is given by (ref. 15) 

A = 4ltR p 2 [l + P - (1 - F) 2 / 3 ] (14) 

where the expressions for P and F differ depending upon the nature of the colli- 
sion (peripheral versus central) and the relative sizes of the colliding nuclei. 
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For the case where R T > R p , we have (ref. 15) 


P = 0.125(iiv) 1/2 f-jj 


- 2 


} ~ P 


- 0.125 


0 .5 ( txv) 1/2 (— - 2 1+ 1 

1 H j | V v 


■ft 




and 


= 0.75(1 - v) 1/2 ^ ~ ^ - 0.1 25 [ 3 ( 1 - v) 1/2 - l] ^ 1 ~ 


with 


v = 


r p + 


and 


p = 


+ R 



(15) 

(16) 

(17) 

(18) 

(19) 


Equations (15) and (16) are valid when the collision is peripheral (i.e., the two 
nuclear volumes do not completely overlap) . In this case, the impact parameter b 
is restricted such that 


R 


T 


R_ < b < R + R 
P TP 


( 20 ) 


If the collision is central, then the projectile nucleus volume completely over- 
laps the target nucleus volume (b < Rip - Rp), and all the projectile nucleons are 
abraded. In this case, equations (15) and (16) are replaced by 


(21 ) 


and 


( 22 ) 


and there is no ablation of the projectile, since it was destroyed by the abrasion. 
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For the case where R p > R T and the collision is peripheral, equations (15) and 
(16) become (ref. 10) 


P = 0.1 25(jiv) 


1/2 1 



- 0.1 25(0 .5 ( v/|i) 



^(1/v)(1 


2 1/2 
t* ) 

3 


.][(2 - „)|, 1 1/ ^ 1 - 


(23) 


and 


F = 0.75(1 - v) 


1/2 1 - 


- 0.125 


/ 3 ( 1 - v) 


1/2 


[1 - (1 - 


2 3/2 ,r 

H > ][i. 


(1 - 


2 , 1/2 

J . 



(24) 


where the impact parameter is restricted such that 


R - R < b < R + R 
p T P T 


(25) 


For a central collision (b < R p - R T ) with R p > R T , equations (23) and (24) become 


P = 


i„ - . , 


, - l 

\ v/ 


1/2 


(26) 


and 


F = [l - (1 - H 2 ) 3/2 ] 


'-I 


1/2 


(27) 


For the excess surface area obtained from equation (14) , the excitation energy is 
given by 


E = A • E < 28 ) 

exc s 

where E s , the nuclear surface energy coefficient (refs. 11 and 15) obtained from the 
liquid drop model of the nucleus, is 0.95 MeV/fm 2 
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Ablation Factors 


Depending upon the excitation energy, the excited prefragment may decay by 
emitting one or more nucleons (protons or neutrons), composites (deuterons, tritons, 
^He, or alpha particles), or gamma rays. The probability a-^j for formation of a 
particular final fragment of type i as a result of the deexcitation of a prefrag- 
ment of type j is obtained from the EVAP-4 computer code (ref. 12) by treating the 
prefragment as a compound nucleus with an excitation energy given by equation (28). 
The final fragmentation cross section for production of the type-i isotope is then 
given by 


a (z. , A . ) = Y a. . a (z.,A.) 

F l i .ID abr j j 


(29) 


where cr a k r (Z j ,Aj ) is obtained from equation (6) or equation (11). The elemental 
production cross sections are obtained by summing over all isotope contributions as 


a (z) = l a <z,A> 

P A F 


(30) 


FRAGMENTATION RESULTS 

As an illustrative application of the theory, element production cross sections 
for fragments of calcium (Z = 20) and heavier elements were calculated for ^£>Fe pro- 
jectiles, at 1 .88 GeV/nucleon incident kinetic energy and collided with stationary 
target nuclei of 12 C, 108 Ag, and 208 Pb. These reactions were chosen for analyses 
because of the availability of experimental data for comparison purposes (ref. 13) 
and because relativistic 56 Fe nuc i e i are among the dominant HZE particles of radio- 
biological significance for manned spaceflight (ref. 2). 

Tables I through III display the elemental production cross sections obtained 
for carbon, silver, and lead targets using both the hypergeometric (eq. (6)) and 
giant dipole resonance (eq. (11)) dispersion expressions. Also displayed are the 
experimental results from Lawrence Berkeley Laboratory (ref. 13). Except for the 
cross sections for Mn production (carbon and silver targets) and V production (silver 
and lead targets), the agreement between theory and experiment is quite good. When 
compared with the predictions obtained with the hypergeometric distribution assump- 
tion of equation (6), the use of the giant dipole resonance expression for charge 
dispersion (eq. (11)) appears to yield slightly improved overall agreement between 
theory and experiment. 

Figures 1 through 3 display the elemental production cross sections obtained 
from equation (11), for the giant dipole resonance dispersion, along with the experi- 
mental data from reference 13. Also displayed, for comparison, are the predictions 
from the semiempirical relations of Silberberg et al. (ref. 14). For the semi- 
empirical relations, the unmodified predictions are displayed. Also displayed are 
the fragmentation cross sections obtained by renormalizing to insure mass and charge 
conservation. Details of the renormalization can be found in reference 20. For the 
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carbon target (ref. 21), y 2 for the giant dipole resonance (GDR) predictions is 
31.6, which is larger than the 19.4 obtained using the Silberberg-Tsao (ST) methods. 
For the GDR, most of the y 2 comes from the Mn overestimate. If that point is 
excluded, y 2 is reduced from 31.6 to 9.4. The comparative results for Ca, Sc, Ti, 
V, and Cr are in better agreement with experiment. For the silver target, the GDR 
X 2 is 57.3 (5.2 if the Mn datum is excluded), whereas the ST y 2 is 32.4 (9.5 if 
the Mn datum is excluded). For the lead target, the GDR y 2 is 4.4 (1.9 if the 
V datum is excluded), compared with the ST y 2 of 52.4 (2.3 if the Mn underestimate 
is excluded). In general, the overall agreement between theory and experiment for 
the abrasion-ablation model is satisfactory considering its simple nature. 

In previous heavy-ion transport work (refs. 22 and 23), it was demonstrated that 
the improved agreement between theory and experiment for Bragg (depth-dose) curves 
was obtained by using Silberberg-Tsao fragmentation parameters modified to scale by 
velocity (rather than total kinetic energy) and renormalized to conserve fragment 
charge and mass. As shown in figures 1 through 3, the modifications (labeled VR) do 
improve the ST predictions for the predominate, near-projectile mass fragments (in 
this case, Mn), but yield substantial overestimates for the fragmentation cross 
sections for the lighter mass fragments. Simple corrections to the Silberberg-Tsao 
parameters, such as renormalization, are apparently adequate for gross total-dose 
comparisons (ref. 23). However, only certain fragments may be biologically signifi- 
cant (ref. 3). Therefore, these corrections may be inadequate for the more pertinent 
shielding problems such as the accurate prediction of individual fragment species 
production. Clearly, the need remains for a comprehensive and accurate HZE particle 
fragmentation theory, of which the work described herein is a beginning. 

Future improvements to this simple abrasion-ablation model will center on 
extending the GDR charge dispersion method to incorporate the actual quantum- 
mechanical abrasion formalism rather than using the geometric model approximation of 
reference 10. Improved methods for estimating the prefragment excitation energy 
spectrum will also be developed. Finally, if necessary, an alternative to the EVAP-4 
ablation code will be tried, such as an intranuclear cascade code (ref. 24); or the 
development of other methods to describe the ablation step will be undertaken. 


CONCLUDING REMARKS 

A simple abrasion-ablation model for describing the nuclear fragmentation of 
galactic heavy ions has been presented. The model was shown to yield reasonable 
agreement with experimental elemental production cross sections. These cross sec- 
tions were obtained at the Lawrence Berkeley Laboratory for the fragmentation of 
relativistic ion nuclei by various targets. Comparisons were also made with predic- 
tions obtained from the semiempirical formulation of Silberberg and co-workers. 

Simple corrections to the latter, such as the mass and charge renormalization which 
previously yielded improved estimates for total radiation dose as a function of depth 
in an extended water column, were shown to be inadequate for accurately predicting 
individual fragmentation cross sections. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 16, 1984 
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SYMBOLS 


A nuclear mass number 

A^ target nuclear mass number 

B(e) average slope parameter of nucleon-nucleon scattering amplitude, fm^ 

b projectile impact parameter, fm 

C(y) Pauli correlation function, defined in equation (5) 

projectile prefragment excitation energy, MeV 

E nuclear surface energy coefficient, MeV 

s 

e two-nucleon kinetic energy in their center of mass frame, GeV 

F function defined in equations (16), (22), (24), and (27) 

GDR giant dipole resonance 

1(b) defined in equation (4) 

J droplet model coefficient (25,76 MeV) 

m number of abraded nucleons 

N total number of nuclear neutrons 

N renormalization coefficient for jth prefragment 

j 

n number of abraded neutrons 

P function defined in equations (15), (21), (23), and (26) 

p(b) probability for not removing a single nucleon by abrasion 

Q droplet model coefficient (11.9 MeV) 

R uniform nuclear radius of projectile, fm 

P 

R uniform nuclear radius of target, fm 

T 

ST predictions of Silberberg-Tsao formulation without renormalization 

u defined in equation (13) 

VR predictions of Silberberg— Tsao formulation with renormalization 

y two-nucleon relative position, fm 

Z total number of nuclear protons 

z number of abraded protons 



longitudinal position of projectile center of mass, fm 
binomial coefficient 

a. . probability of formation of type-i fragment as a result of deexcitation of 

type-j fragment 

a defined in equation (12) 

z 

p defined in equation (18) 

A excess nuclear surface area, fm 2 

\x defined in equation (19) 

v defined in equation (17) 

collection of constituent relative coordinates for target, fm 
p nuclear density, fm” 3 

p p projectile nuclear density, fm” 3 

P T target nuclear density, fm" 3 

a(e) average nucleon-nucleon total cross section, mb 

a abr (Z,A) cross section for production of nucleus of type (Z,A) by abrasion, mb 
fragmentation cross section, mb 
cross section for abrading m nucleons, mb 
An arrow over a symbol denotes a vector* 
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TABLE I.- ELEMENTAL PRODUCTION CROSS SECTIONS FOR THE REACTION 

56 Fe + 12 C -► Z + X 

[Incident kinetic energy is 1.88 GeV/nucleon] 


Elemental production cross sections, mb 


Element 

produced 

Hypergeometric 

Giant dipole 
resonance 

Experiment ( ref • 13) 

Fe 

161 

209 


Mn 

321 

308 

181 ±27 

Cr | 

156 

142 

124 ± 13 

V 

126 

124 

100 ± 11 

Ti 

90 i 

88 

87 ± 11 

Sc 

69 

69 

54 ± 9 

Ca 

77 

78 

78 ± 11 


TABLE II.- 


ELEMENTAL PRODUCTION CROSS SECTIONS FOR THE REACTION 
56 pe + 108 Ag + z + x 


[Incident kinetic energy is 1 .88 GeV/nucleon] 


Element 

Elemental 

production cross sections, mb 




produced 


Giant dipole 



Hype r ge ome trie 

resonance 

Experiment (ref. 13) 

Fe 

296 

262 


Mn 

381 

446 

280 ± 23 

Cr 

226 

230 

218 ± 21 

V 

150 

149 

117 ± 15 

Ti 

126 

128 

124 ± 16 

Sc 

101 

100 

104 ± 13 

Ca 

j 

102 

112 

118 ± 14 


TABLE III 


ELEMENTAL PRODUCTION CROSS SECTIONS FOR THE REACTION 
56 Fe + 208 Pb -► Z + X 


[Incident kinetic energy is 1.88 GeV/nucleon] 


Element 

produced 


Elemental production cross 


Hypergeometric 


Fe 

345 

Mn 

445 

Cr 

267 

V 

175 

Ti 

152 

Sc 

121 

Ca 

116 


Giant dipole 
resonance 

302 

521 

268 

174 

151 

119 

129 


sections, mb 


Experiment (ref. 


509 ± 40 
242 ± 25 
142 ± 20 
148 ± 22 
111 ± 17 
144 ± 22 


13) 


Production cross section, mb 



Ca Sc Ti V Cr Mn Fe 

Element produced 


Figure 1 •- Elemental production cross sections for iron projectile nuclei 

fragmenting in carbon targets. 
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production crobb bccuiun 



Figure 2.- Elemental production cross sections for iron projectile nuclei 

fragmenting in silver targets. 



Production cross section, mb 



Element produced 


Figure 3.- Elemental production cross sections for iron projectile nuclei 

fragmenting in lead targets. 
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